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Abstract 



We study the spontaneous symmetry breaking of 0(3) scalar field on a fuzzy sphere 
Sp. We find that the fluctuations in the background of topological configurations are 
_ , finite. This is in contrast to the fluctuations around a uniform configuration which 

Qj _ diverge, due to Mermin-Wagner-Hohenberg-Coleman theorem, leading to the decay 

j-S^ ■ of the condensate. Interesting implications of enhanced topological stability of the 

configurations are pointed out. 
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1 Introduction 

cn_ 

^ ■ Field theories on non-commutative geometries are inherently non-local[Tl El [3l HI [5] . This 

feature gives rise to novel behaviour such as the mixing of infrared and ultraviolet scales. 
These lead to new ground states with spatially varying condensates |5l 13 El IH]- Many 
non-perturbative studies have established that fuzzy spaces such as Groenewold-Moyal 
plane, fuzzy spheres allow for the formation of stable non-uniform condensates as ground 
[Xj ■ states[IOl [niliainillllllSllISllIIlIIHllIH]. NonlocaUty plays an essential role in all these. 

Exploring implications of this nonlocal nature of field theories is one of the intensive research 
activity in the field of noncommutative physics today. 

In the past, we have analysed the spontaneous breakdown of global symmetries and 
the corresponding behaviour of Goldstone modes through simulations |17[ [T8] . The main 
aim of this paper is to study the interplay of SSB and topological features of the sigma 
models on fuzzy spheres. This issue is important because the Mermin-Wagner-Hohenberg- 
Coleman (MWHC) theorem states that there can be no SSB of continuous symmetry on 
2-dimensional commutative spaces [20 | UT ] 12 ^ [25]. This theorem was established rigorously 
by summing up all the infrared contributions in the context of field theories with local 
interactions. The resulting diverging fluctuations of the Goldstone modes prevent the sym- 
metry from being spontaneously broken. No SSB also implies no magnetisation. There is 
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no obvious generalisation of the MWHC theorem for the non-commutative spaces, since the 
theorem rehes strongly on the locality of interactions. Gubser and Sondhi, [B], have shown 
that the Goldstone mode fluctuations become severe in the non-commutative case. These 
authors infer, after analysing the fluctuations in the background of a uniform condensate, 
that a condensate with only zero momentum/ angular momentum mode are unstable. How- 
ever non- commutative spaces admit non-uniform solutions and one can ask the question 
what happens to the stability of these configurations. Non-uniform condensates naturally 
have an infra-red cut-off for the fluctuations. This cut-off can soften the otherwise divergent 
contributions of the Goldstone modes. It would be desirable to compute the fluctuations by 
analytical methods as was done for the MWHC theorem. Only few models allow for exact 
treatment of SSB even in the commutative spacetime. It is difficult to analytically sum the 
fluctuations in the present case. But in the absence of such a study numerical simulations 
can provide answers in the present context. 

In a previous work of ours, we explored these issues for U{1) symmetry using numerical 
simulations [18]. As expected, we found that the uniform condensates are unstable. But 
the condensates with higher angular momentum components still survived the fluctuations. 
However when the thermodynamic and continuum limits were taken, the flutuations in 
the momentum mode above the cut-off increased. The Goldstone modes were capable of 
destroying the low lying modes of the non-uniform condensate. 

In this paper we analyse the 0(3) symmetry on fuzzy sphere. Nonlinear sigma models 
on fuzzy spheres have been studied in the literature exploring topological features |2H ES] . 
We expect that if we repeat the calculations of U{1) for this case the results will be more 
or less similar. But the 0(3) case is significantly different from our previous study in an 
important way. In this case there are non-uniform configurations which are topologically 
stable. For example the hedgehog type configuration. A condensate which is topologically 
stable can be unstable against fluctuations with energies of the order of the system size. So 
we expect that Goldstone modes whose energy goes as logarithm of the system size may not 
be able to destroy the hedgehog condensates. To settle this, we calculate the fluctuations 
around such configurations via numerical simulations. Our numerical results show, that 
indeed, topological stability and nonlocal interactions make the lowest non-zero mode of 
the condensate stable. 

This paper is organised as follows: In Sec. 2 we describe the model and the non-trivial 
topological configurations. The numerical simulations and the results are discussed in Sec. 3. 
In Sec. 4 we provide our conclusions. 

2 0(3) Model and Topological Condensates on 5*1^ 

The fuzzy sphere is described by the coordinates Xj satisfying the following SU (2) algebra, 

N 

Here R is the radius of the fuzzy sphere. One can choose the angular momentum j ~ ~^ 
representation for the fuzzy sphere. 
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With this choice the most general action for 0(3) hermitian scalar field $j (i = 1, 2, 3) 
upto quartic interactions takes the form, 



Air 

^($) = -Tr 



(2) 



In the mean-field the above theory admits a uniform condensate for r < 0. However 
as discussed above the fiuctuations of the Goldstone mode render this solution unstable. 
Apart from the uniform condensate the above model admits many meta-stable solutions. 
To simplify our arguments we first consider the case P = 0,fi = 0: 

The number of meta stable configurations increases with the matrix size of $. Among 
them we are interested on those which are stable due to topological obstructions. For 
example, 

$i = aLi, with a = y ^ (3) 

The analog of this configuration in continuum space is the hedgehog configuration where 
the 0(3) spin vector on the sphere is pointing radially outward. The spin is parallel to 
the position vector on the sphere. This configuration is topologically stable as it cannot be 
smoothly deformed to a uniform one. Similarly the above configuration cannot be smoothly 
deformed to $ = / which is also a solution. The above configuration corresponds to a 
winding one map from the physical space Sp to the vacuum manifold which is Sp. All 
topologically stable configurations in the continuum limit, can be characterised by the 
second homotopy group 112(5'^). For a discussion on topological classification of the maps 
Si see [21 Eg. 

When we include the other terms in the action we have more parametric choices for the 
topological configurations. For example when (3^0, but /i = we have 



$i = aLi, with a = ^ ^ 

Then Eq.(jlj) shows clearly the possible existance of topological vacuum configurations 
even when r > 0. For example, when r = we have 

3/3 

= aLi, with a = ^ _ (5) 

While it is important to study the topological stability in such a general situation the results 
are not characteristically different from /3 = 0. Here we focus on /3 = in this paper for 
simulations but discuss general cases. 

To study the net effect of topological nature of the background configuration and non- 
locality on fiuctuations, we consider only the winding number one configuration which is 
given in Eq.([3]). As mentioned our plan is to compute these fiuctuations numerically. Even 
before computing the fluctuations one can make some general remarks about the behaviour 
of the fluctuations [18]. The effect of nonlocality basically provides a non-zero mass ^ij^) 

to the Goldstone mode fluctuations. This puts an infrared cut-off for the fluctutations. 
From our previous study [18] it seems that this mass/cutoff is mode dependent, as only 
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higher modes of the condensate survived the fluctuations. As we will see from our results 
the combined effect of topology and nonlocality, the infrared cut-off drastically reduce the 
contribution of the fluctuations. In the next section we describe our numerical simulations. 



3 Numerical simulations 

We use "pseudo-heat bath" updating method in our numerical simulations which is de- 
scribed in detail in our earliar papers [171 IE] ■ In our method, for each choice of parameters, 
we choose an initial conflguration given by Eq. (|3]). Fluctuations around this conflguration 
are then generated by the above updating method. Since this conflguration is a variational 
solution to minimising the classical action, it will thermalise as we update/include the ther- 
mal fluctuations. Once the initial conflguration is thermalised we compute the observable 
M. We make measurements after every 10 updates of the entire matrix. We also use over- 
relaxation to reduce the auto-correlation of the conflgurations generated in the Monte-Carlo 
history. 

In a numerical simulation, the condensate will not maintain its exact form as in Eq. 
along the Monte-Carlo history. The conflguration can evolve into different random SU (2) 
rotated conflgurations of Eq.Q as we keep updating it. To overcome this, one needs to 
rotate the conflguration at each step of the Monte-Carlo history so that the conflguration 
takes the form of Eq. (|3]). But this is a difficult and time consuming task. On the other 
hand one can have an observable made of $^s which is invariant under the SU{2) rotations, 
e.g basis independent. For this purpose we define the following observable. 



A,, = — Tr(L,$,),M = v^ (6) 

M projects out the I = 1 angular momentum mode. Note that the initial configuration in 
Eq. ([3]) projects out only the I = 1 mode. Analysing the statistical behavior of M will 
give us a definite conclusion about the stability of the initial configuration. We mention 
here that Tr(^^ is also an SU{2) invariant. But the information on the amplitudes of 

i 

different / modes gets lost in this form. Also comparatively the observable M may serve 
as an order parameter in the case of any phase transition of the hedgehog configuration 
to $i = at high temperatures. We mention here that / = 1 is the lowest possible stable 
mode, as / = mode will be unstable. One can consider configurations with higher winding, 
instead of Eq.([n]), however we expect them to be more stable than the / = 1 condensate. 
This is because the infrared cut off will rise with higher winding configurations. 

For practical reasons, the size of the matrix A^, in other words size of the resolution scale 
is finite. So there are usually finite volume (i?, A^) effects. So a non- vanishing condensate 
$j does not mean there is SSB. One needs to define suitable observable dependent on $j 
which should scale with A^) appropriately in the thermodynamic limit (A^ — )■ oo, i? — > 
oo) to conclude anything. Now there are two possible thermodynamic limits. If in the 

thermodynamic limit the ratio — does not vanish then the space is described by a non- 
commutative algebra. This limit is of interest to us, as we expect that the MWHC theorem 
will hold good in the commutative thermodynamic limit. 
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Figure 1: M vs Monte Carlo history for = 48 
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Figure 2: Histogram H(M) for = 48 

4 Results and discussions 

In our calculations we fix — = 10, r = —8. For simplicity we take Ai = 0.25. With 

this choice of parameters we do our simulations for five different sizes of the $ matrices, 
= 48,56,64,78,96. FIG.l gives a typical Monte-Carlo history of our simulations for 
= 48. 

In FIG.l M fluctuates around a value close to the initial value. Then M suddenly 
jumps to a small value and settles down. A histogram H{M) of M clearly shows two peaks 
M as seen in FIG. 2. The peak on the left has large / = and small / = 1 component. 
The peak at higher value of M has large / = 1 component and small / = 0. This peak 
is close to the value of the initial configuration. So in this state fluctuations modify the 
initial configuration slightly and retain its topological nature. In our Monte-Carlo history 
we observed the I = 1 state decaying to / = state but not vice-versa. This implies that 
due to finite volume effects, the uniform condensate is more stable than our initial hedgehog 
configuration for this case of A^ = 48. 

To study the stability of the / = 1 configuration we considered both the commutative 
and non- commutative limit. For the commutative limit we fixed and considered higher 
values of A^. We did not observe any change in the distribution of M in the / = 1 state. 




H(M) 
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The average value, and the fluctuations of M remain almost the same as we go from 
= 48 — )■ 64, as can be seen in FIG. 3. As for = 48 the / = 1 configuration also decays 
for A^ = 64. This result suggests that the / = 1 topological configuration is not stable in 
the commutative continuum limit as expected. 

1.2 I , , , , 1 
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Figure 3: M-C history for fixed 
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Figure 4: M-C history for fixed R^/N 

There is a complete change in the behavior as we consider the non-commutative limit, 

i.e fixed — as we increase A^. Except for the lowest A^ = 48 the I = 1 state did not 
A^ 

decay during the entire run for higher A^. In FIG. 4 we show the Monte-Carlo history of 
A^ = 48, 64, 96. Unlike the commutative limit, the fluctuations of M decrease with A^. This 
can be clearly seen in FIG. 4. In FIG. 5, we give the average value of M as a function of 
A^. The average value of M increases slightly with A^, with the variation decreasing with 
A^. This suggests M will reach a finite value in the continuum limit. We also compute 
the fluctuations of M to see any possible scaling with the cut-off A^. In FIG. 6, we show 
X = (M^) — (M)^ in the / = 1 state. The solid curve represents a fit, f{N) ~ A^" with 
a ~ —4.. This clearly suggest that the I = 1 state is stable in the N ^ oo leading to 
spontaneous breaking of the 0(3) symmetry. 

We also mention here that one can start with an initial uniform / = configuration and 
consider fluctuations. We expect that the results be similar to that in ref.[T8]. In ref.flH] 
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it was found that only the highest mode / = — - — condenses. The fact that we find 

the / = 1 mode stable clearly shows that the topological nature of the initial configuration 
complements the effect of non-locality. These two effects drastically reduce the fluctuations. 



5 Conclusions 

In this analysis we have shown that topologically non-trivial configurations on the fuzzy 
sphere avoid the MWHC theorem much more dramtically than the non-topological sym- 
metry breaking. The mass gap or the infrared cut-off in this case is large enough to render 
the fluctuations of the Goldstone modes finite. On the other hand for non-topological con- 
densates the Goldstone modes are large enough to destroy almost all the modes except the 
few highest modes. It is interesting to scale A and r as we go to continuum limit. The limit 
of noncommutative continuum geometry is always expected to maintain SSB and stable 
solitons as there is obstruction to MWHC theorem from nonlocality. This analysis will be 
presented elsewhere. 

We have presented the simulations wherein the cubic Chern-Simons (CS) term is absent 
in the action Eq ([2]). The Chern-Simons term allows topological solitons even when the 
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quadratic mass term is positive upto some value. Interestingly with CS term the config- 
uration (pi = a Li is preferred over symmetric solution. On the otherhand, it is not 
expected to alter the picture of topological stability of the solutions. This term plays an 
important role in the emergent geometry in NC fuzzy spaces [26l |27] . What we find here 
in the simulations is that even in the absence of CS term, emergent fuzzy spaces can be 
stable. The stabilty of higher dimensional fuzzy spaces like CPp are of significance in this 
context [2B]- The implications of this stability for extra-dimensional fuzzy spaces will be 
considered later. 
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